Long-range interacting many-body systems with alkaline-earth-metal atoms 
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Alkaline-earth-metal atoms exhibit long-range dipolar interactions, which are generated via the 
coherent exchange of photons on the 3 Po— 3 Di-transition of the triplet manifold. In case of bosonic 
strontium, which we discuss here, this transition has a wavelength of 2.7 fim and a dipole moment of 
2.46 Debye, and there exists a magic wavelength permitting the creation of optical lattices that are 
identical for the states 3 Po and 3 Di. This interaction enables the realization and study of mixtures 
of hard-core lattice bosons featuring long-range hopping, with tuneable disorder and anisotropy. We 
derive the many-body Master equation, investigate the dynamics of excitation transport and analyze 
spectroscopic signatures stemming from coherent long-range interactions and collective dissipation. 
Our results show that lattice gases of alkaline-earth-metal atoms permit the creation of long-lived 
collective atomic states and constitute a simple and versatile platform for the exploration of many- 
body systems with long-range interactions. As such, they represent an alternative to current related 
efforts employing Rydberg gases, atoms with large magnetic moment, or polar molecules. 



Introduction.- The creation and exploration of quan- 
tum systems with long-range interactions is in the focus 
of intense research activity worldwide. Within the con- 
text of novel technological applications, such as quantum 
information processing, strong long-range interactions 
are essential as they permit the implementation of en- 
tangling gate operations among distant qubits. From the 
perspective of fundamental physics of condensed matter 
systems, these interactions permit the study of strongly 
correlated phases of quantum matter. In order to access 
this potential there is a need for a simple experimen- 
tal platform that fosters long-range interactions. In the 
domain of ultra cold gases there are currently three ap- 
proaches, which rely on atoms with large magnetic dipole 
moment (e.g. chromium, dysprosium and erbium PUS]), 
polar molecules [4] or atoms excited to Rydberg states. 
In particular Rydberg atoms have celebrated recent suc- 
cesses in the realms of fundamental physics and techno- 
logical applications [SHE] • Two experiments have success- 
fully implemented quantum gate protocols among qubits 
encoded in distant atoms [9] [10] and, very recently, the 
intricate dynamics of strongly correlated Rydberg lattice 
gases was studied in experiment [TTHT3] . 

In this work we describe a novel platform for the real- 
ization of many-body systems featuring long-range dipo- 
lar interactions. It is based on the exchange of virtual 
photons between low-lying triplet states of alkaline-earth- 
metal atoms. We specifically focus on bosonic strontium 
(Sr) atoms trapped in a deep optical lattice in a Mott 
insulator state and photons exchanged on the 3 Po— 3 Di- 
transition (see Fig. [IJi). We derive the corresponding 
many-body Master equation, focussing specifically on the 
situation of planar and linear x^-models with long-range 
interactions, which are equivalent to hard-core lattice 




FIG. 1. a: Relevant levels of the Sr atom, b: Atoms are 
trapped in an optical lattice and the interaction between them 
is generated by the exchange of (virtual) photons on the tran- 
sition between 3 Po and the three degenerate 3 Di states. The 
decay rate and wavelength are T = 290 kHz and A = 2.7 /mi, 
respectively. A is much larger than the typical interatomic 
spacing (a = 206.4 nm) at the magic wavelength. 



bosons with long-range hopping. We provide data of a 
magic wavelength for an optical lattice that grants equal 
confinement for both the states 3 Po and 3 Di, characterize 
the role of decoherence and disorder and show how the 
interactions and the resulting collective light-scattering 
become manifest in the fluorescence spectrum. Build- 
ing on the routine creation of alkaline-earth-metal Mott 
insulators in a number of laboratories [14] [15] , our ap- 
proach represents a simple route for the exploration and 
exploitation of many-body phenomena in long-range in- 
teracting systems and highlights a novel way for the cre- 
ation of long-lived collective atomic states, with applica- 
tions in quantum optics and quantum information. 

Triplet states of strontium.- Sr has two valence elec- 
trons and its spectrum is thus formed by a series of sin- 
glet and triplet states (Fig. [Ik). The radiative transitions 
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between the two series are dipole forbidden, which - due 
to the resulting small transition line widths - leads to a 
wide range of applications, such as ultra precise atomic 
clocks [16] or the implementation of quantum informa- 
tion processing protocols [5j [17]. Here, we consider Sr 
atoms in a Mott insulator state [14] [15] as depicted in 
Fig. [TJd. The lattice is identical for the 3 P and 3 Di 
states, and its blue-detuned magic wavelength is located 
at Afr m = 412.8 nm (for more detailed information, see 
Supplementary Material). The resulting lattice spacing 
is a = 206.4 nm, a value that we will use throughout this 
work to benchmark our results. 

Initially all Sr atoms within the lattice are excited to 
the triplet manifold, i.e. either to the metastable state 
3 Po (its lifetime can be considered infinite for all exper- 
imental purposes) or the state 3 Di. Long-range interac- 
tions between two Sr atoms then emerge by the resonant 
exchange of photons that are emitted on the transition 
3 Di— 3 P which has the wavelength A = 2.7 pm (see Fig. 
[TJd). This mechanism, which leads to a dipolar interac- 
tion, is in general well-understood [T8H2Q] but usually it 
plays a role only in very dense samples, i.e. where the 
interatomic distance is far smaller than the wavelength of 
the transition [21] . Conventional lattice setups usually do 
not reach such parameters. However, the combination of 
a lattice with small spacing and a long wavelength transi- 
tion that is available in Sr allows us to enter this regime. 
The transition dipole moment between the 3 Po state and 
the three degenerate 3 Di-states is p = 2.46 Debye and ef- 
fectuates a strong resonant dipole-dipole interaction that 
extends over several lattice sites as depicted in Fig. ^p. 
In order to avoid unnecessary complications we restrict 
ourselves to photon exchange on the 3 Di— 3 Po transition, 
which has the highest decay rate (r = 290 kHz, [22] ) 
and leave the (straightforward) consideration of the addi- 
tional weaker coupling channels 3 Di— 3 Pi and 3 Di— 3 P2 
to future investigations. 

Many-body Master equation.- The starting point for 
the derivation of the many-body Master equation is the 
Hamiltonian describing N Sr atoms coupled to the ra- 
diation field. To formulate it we introduce the vec- 
tor transition operator for the k-th atom (located at 
r /e) bfc = b kx x + b ky y + bkz z such that the transi- 
tion dipole matrix elements are real and aligned along 
the three cartesian spatial axes x,y and z [23]. Here, 
bkj = \P)k 01 w ^ n j — %,y,z, where \P) k represents the 
k-th. atom in the 3 Po state and \j) k the cartesian states 
of 3 Di, related to the angular momentum ones \m) k 
(with m = -1, 0, +1) as |=pl) = (± \x) - i \y)) /y/2 and 
|0) = \z). Within this notation the Hamiltonian of the 
atomic ensemble and the radiation field is given by H a ? = 

EfcLl ^ ab ! ' hk + ^qA ^q«qA a qA + EfcLl E q A gqA ' 

( a q A s fc e " iq ' rfc - s/ca qA e iq,rfe ) with s k = b[ + b fc , where 

huj a = 2ithc/\ is the energy difference between the 3 Po 
state and the degenerate 3 Di-manifold, fkj q is the en- 



ergy of a photon with momentum q and polarization A 
and a q A is the annihilation operator of such a photon 
(bosonic operators, i.e., [a q A,a q , A ,] = Sq^Sw). The co- 
efficient g qA is given by g qA 



P\j2^y^ with V the 
quantization volume, and e q A the unit polarization vector 
of the photon (q • e q A = 0). 

Following Refs. [20 , 23 we obtain the Master equation 
governing the evolution of the density matrix p of the 
atomic ensemble: p = — | [H,p] -\-V(p). The first term 
depends on the many-body Hamiltonian 

H = hu; a bl ' h k + ft b l ' Wi ' hi. (1) 

k k^l 

Its first part contains the bare energies of the atomic lev- 
els and the second part describes the long-ranged and (in 
general) anisotropic dipole-dipole interaction, character- 
ized by the coefficient matrix 



Vki 



yo (km) 



yi (km) 



km 



Here y n (x) represents the n-th order spherical Bessel 
function of the second kind and n k i = k a rki with fc a = 
uj a /c and r^i = — r/ = VkiVki- The second term of the 
Master equation depends on the dissipator 
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J2bk-r kl .pb}-Ub{.r kl .b hP }. 

kl 



The coefficient matrix 



kl 



jo (km) 



ji (km) 



km 



encodes the dissipative couplings among the atoms and 
j n (x) represents the n-th order spherical Bessel function 
of the first kind. 

The coherent and dissipative dynamics are intimately 
connected since they both originate from the emis- 
sion/absorption of photons. However, by virtue of the 
long wavelength of the photons and the achievable small 
interparticle separation (see Fig. [TJd) one can reach a pa- 
rameter regime in which the coherent interaction is much 
stronger than the dissipation. This is shown in Fig. [2^i 
where we compare the coherent interaction V\2 to the dis- 
sipative rate Y\2 for two atoms separated by a distance 
d whose induced dipoles are aligned as tt- For d = a 
the ratio Vyil^vi is approximately 6.5, but for the dipole 
alignment this ratio reaches 15.4, showing that the 
Sr lattice setup is indeed well-suited for the study of co- 
herent many-body phenomena. In fact, we will later see 
that due to the presence of sub-radiant states even larger 
ratios can be achieved. 

Coherent dynamics.- Hamiltonian conserves the 
number of atoms in the 3 Di-state, TVd = J2k ^1 ' 
Hence, non- trivial dynamics is only induced by the sec- 
ond term of eq. ([!]) which depends strongly on the geom- 
etry of the lattice. To get a glimpse of the versatility of 
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FIG. 2. Transport of a single excitation on a chain of N = 20 
atoms with induced dipoles pointing in the ^-direction, a: 
Comparison between the nearest neighbor interaction V12, the 
damping rate T12, and the effective decay rate r e ff of the 
single excitation for different values of the nearest neighbor 
separation d. For d — a the effective decay rate is approx- 
imately two orders of magnitude smaller than the coherent 
interaction, b, c, d: Time-evolution of a single, initially lo- 
calized excitation, for varying disorder and for d — a with: b: 
a/a = 0, c: a/a = 0.025 and d: a/a = 0.05. 



the Sr setup for the study of coherent many-body phe- 
nomena let us consider a situation in which we have a 
one- or two-dimensional lattice located in the x — y plane 
and where atoms are solely excited to the 3 Di(m = 0)- 
state. In this case the Hamiltonian ([!]) simplifies to 
H xy = ^2 k ^i Wkib\ z bi z with the interaction coefficients 
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yo (km) 



yi(Kki) 
km 



(the dipo- 



4fef ^ |r fc -n| 3 

lar approximation holds when n^i < 1). H xy represents 
an xy-m.ode\ with long-range interaction or equivalently 
a system of hard-core bosons with long-range hopping, 
which is due to the fact that the operators b\ z and bk z 
can be identified with creation/annihilation operators of 
hard-core bosons. Note, that such long-range hopping 
can also be effectively established in ion traps [24] [25] . 
For nearest neighbor interactions the xy-mode\ has been 
studied extensively in the literature in the context of 
quantum information processing [26] , quantum and ther- 
mal phase transitions [27] and relaxation of closed quan- 
tum systems [28 . The case of spin-systems with long- 
range interactions is less explored [29] and recent Monte 
Carlo simulations [30] raise new questions concerning 
their phase behavior. 

To study the (thermo) dynamics of H xy one needs to 
control the density of hard-core bosons N-&/N which is 
done as follows: Starting with all atoms in the state 
3 Po one irradiates a laser on the 3 Po- 3 Di-transition. 
The Hamiltonian describing the atom-laser coupling is 

h l = nn L • Y^k=i (e _i(k ' rfc_u;£) b fe + e *(k-r fe -u;t) b t^ with 

= pEo/ft, where uj is the frequency, k the momen- 
tum and Eq the amplitude vector of the laser. Apply- 



ing a strong laser pulse | ^> with amplitude 

vector E = Eqz for a time r on resonance (u = u a ) 
creates on average Nn(r) = N sin 2 (2\fti,\r) hard core- 
bosons. This number fluctuates with standard deviation 
A7V D = v / iV|sin(4|fi L |r)|/2. Alternatively, one could 
think of selectively changing the state of atoms in spe- 
cific sites. In particular, in the one-dimensional case, the 
required single-site addressing can be achieved by apply- 
ing a magnetic field gradient, which is switched off after 
the desired state is prepared [31] [32] . 

The xy-mode\ focussed on here merely represents a 
very simple scenario. In all its generality the Hamil- 
tonian ([I]) describes three species of hard-core bosons 
which, depending on the lattice dimension and geome- 
try, exhibit anisotropic and species-dependent long-range 
hopping and species inter conversion. The density of the 
individual species is furthermore controllable by simple 
laser pulses. This generates a rich playground for the 
discovery and study of many-body quantum phases. 

Dissipation and disorder.- We now discuss two seem- 
ingly harmful effects: Collective dissipation due to radia- 
tive decay and disorder stemming from the uncertainty 
in the atomic positions. To analyze them we numerically 
simulate a situation in which a single excitation (or hard- 
core boson) propagates through a linear chain of N = 20 
Sr atoms orientated along the x-direction under the ac- 
tion of H xy and the corresponding dissipator. Initially, 
the excitation is localized at the leftmost atom which is in 
the 3 Di(m = 0)-state. Its propagation is depicted in Fig. 

which also shows an expected overall decrease of the 
excitation density due to dissipation. Moreover, we have 
extracted from the simulations the (effective) decay rate 
r e ff of the excitation as a function of the lattice spacing 
d (see Fig. [2|l). Surprisingly, one can see that for d — a 
the effective decay rate of the excitation is not only much 
smaller than the single-atom one 7 = T/3, but more than 
two orders of magnitude smaller than the dipole-dipole 
interaction. The reason for this unexpected long lifetime 
of the excitation resides in the collective character of the 
dissipation and the emergence of sub-radiant states. 

Let us now discuss disorder, which arises from the 
fact that the external wavefunction of each localized 
Sr atom has a finite width. We model this wave- 
function as a three-dimensional isotropic Gaussian with 
width o~ which is localized on the respective lattice 
site. Using the dipolar approximation for the inter- 
action and assuming r k i ^> o~ the couplings Wki be- 
come random variables which are distributed accord- 
ing to p(W k i) = [(A 2 / 3 )/(3V2^cir kl W^ /3 }e W [-(rki ~ 
[A/W k i] 1/3 ) 2 /(2a 2 )], where A = KT/(4kl). Again, the 
simulated excitation transport nicely reveals the effect of 
disorder in Figs. [2Jd, c and d for d = a. With increasing 
ratio a /a transport becomes less efficient and more popu- 
lation remains localized at the left of the lattice. Beyond 
this very simple illustration it will be very interesting to 
study the effect of this controllable disorder in the many- 
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FIG. 3. a: One-dimensional lattice (oriented along the x axis) 
with lattice constant d. A laser is shone on the atoms with 
momentum k and polarization Eo both perpendicular to the 
chain, b: Power spectrum of the radiation scattered into the 
y — z plane S(u) (arbitrary units and logarithmic scale) as a 
function of the interparticle separation d and the frequency of 
the laser uj. The two insets show the corresponding spectrum 
for d = a and d — 2a (linear scale). The average excitation 
number on the j-th lattice site Nr>j, corresponding to the 
states represented by two of the peaks is also sketched. 



body context. It has been shown that the ground state 
of hard-core bosons exhibits a localization transition for 
certain types of disorder [33] in the hopping rates. It is an 
open question whether this transition is actually present 
also here. 

Spectroscopy on the few-body system.- Finally, let us 
discuss the spectroscopic properties of the Sr lattice 
system and find out whether they provide a clear ex- 
perimental signature of the presence of long-range in- 
teractions. To this end we calculate the power spec- 
trum of the radiation field E(R, t) scattered off the 
ensemble of atoms in the direction of observation R 
driven by a weak incident laser field of frequency 
to: S(R,u>) = £/^dre <WT <Et(R,t).E(R ) t + r)). 
We consider again the above-discussed case of a one- 
dimensional lattice with N sites oriented along the x-axis 
where nearest neighbors are separated by the distance d. 
The dynamics of this driven ensemble is described by 
Hamiltonian H xy , the corresponding dissipator and the 
Hamiltonian i^L, which takes into account the action of 
the laser field with polarization Eo = Eqz. For the sake 
of simplicity we orient the laser beam such that its mo- 
mentum k is perpendicular to the chain (see sketch in 
Fig. [3^t) and observe the spectrum of the radiation scat- 
tered into the y — z plane (denoted by S(u)). 

Before discussing the results for a chain, let us first con- 
sider the case of two atoms, i.e. N = 2, where the prob- 
lem can be treated analytically. Here, the spectrum of 
the radiation S(uj) consists of a single peak of Lorentzian 
form S(lj) oc 1/[(cj - w a - W 12 /K) 2 + (T/6 + T 12 ) 2 ] with 
Ti2 = j\jo (kd) — ] • Hence, the signature of the 
strong interaction is a shift towards the blue (repulsive in- 
teraction) and broadening of the Lorentzian peak [34| [35] . 



Note that only one peak (the symmetric superposition of 
the two singly-excited states) appears in this case. This 
is due to two facts: (i) In this particular geometry (the 
momentum of the laser being perpendicular to the inter- 
particle separation) the laser only couples to the (super- 
radiant) symmetric state and (ii) the symmetric and anti- 
symmetric states are eigenstates of the dissipator as well 
as of the Hamiltonian and, hence, the dissipation induces 
no couplings between the symmetric and anti-symmetric 
states. 

For a chain of N = 20 atoms we have calculated the 
spectrum numerically and the result is shown in Fig. (3}d. 
One observes that, as in the case of two atoms, the laser 
field couples to the symmetric (spin wave) state, which 
decays super-radiantly (visible as the very broad feature 
on the blue side of the atomic line). More interesting, 
however, is the emergence of a number of narrow peaks 
belonging to long-lived sub-radiant states [36]. These 
lines appear due to the fact that the coherent and dis- 
sipative interaction do not share the same set of eigen- 
states. This admixes a fraction of the symmetric state 
to other collective states, making them excitable by the 
laser. The density profile of two selected states is shown 
in the inset of Fig. [3]o [37 . The Sr system thus offers the 
possibility of exciting long-lived collective states [38] [39] 
which can find an application in quantum information 
and photon storage. 

Conclusions and outlook.- In conclusion, we have 
shown that alkaline-earth-metal atoms in optical lattices 
provide a platform offering controllable many-body sys- 
tems with long-range interactions. Specifically, we have 
analyzed a regime which implements hard-core bosons 
with coherent long-range hopping. Beyond that simple 
example, the system permits the study of mixtures of 
hard-core bosons in the presence of tunable disorder. Sig- 
natures of the long-range interaction are manifest in the 
spectrum of the radiation which is collectively scattered 
from the atomic ensemble. In the future we will investi- 
gate closer the properties of the emitted light [40] and its 
use to detect quantum phases [41 ] and (disorder-driven) 
phase transitions. We will furthermore investigate dy- 
namical phases [42] resulting from the interplay between 
dissipative and coherent dynamics and explore how the 
versatility of the Sr platform can be further enhanced, 
e.g., by the application of static electric and microwave 
fields. 
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Supplementary Material for "Long-range interacting 
many-body systems with alkaline-earth-metal 
atoms" 



In this supplementary material, we provide some de- 
tails on a possible implementation of an optical lattice 
configuration for strontium atoms, which grants equal 
confinement for the (5s5p) 3 Po state and all magnetic 
sub- levels of the (5s4i) 3 Di state, thereby realizing a low- 
dimensional long-range interacting many-body system. 

We consider a one-dimensional (ID) blue-detuned op- 
tical lattice along the x-axis. The positive frequency 
part of the lattice field is given by Ei^~\x)e~ ZUJbt = 

£ b (y[ly + \J\%) cos (^fj e~ iuJbt . Here, £ b is the laser 
field amplitude, uu^ = 2i\cjX\ ) is the laser frequency and c 
is the speed of light in vacuum. Following Ref. [43 , one 
finds that with this choice of the lattice field all three 
states (5s4cf) 3 Di(m = 1,0,-1) have the same ac polar- 
izability a b D (X b ). 

Based on the data for wavelengths and Einstein coef- 
ficients for the electric dipole transitions relevant to the 
(5s5p) 3 P and (5s4<f) 3 Di states in Refs. [44H47] . we have 
calculated a b D (Xi)) and also the ac polarizability ap(A&) 
of the (5s5p) 3 Po state as a function of the lattice wave- 
length In Fig. [Slj^a) we show this data in the vicinity 
of 400 nm. A magic wavelength is found at Afo m « 412.8 
nm where a b P (X} )rn ) = a b D (Xb rn ) « —0.76 x 10 3 a.u. 

Note that the negative polarizability means that 
atoms will be trapped in the intensity minima of the 
lattice field. Thus, in order to tightly trap atoms 
in the ID blue-detuned optical lattice in experiment, 
one needs to provide additional transverse confine- 
ment. Without introducing a differential light-shift be- 
tween the 3 Po and 3 Di states, this additional con- 
finement can easily be created with a two-dimensional 
(2D) red-detuned optical lattice near a red magical 
wavelength. The positive frequency part of the cor- 
responding laser field is given by J^~\y,z)e~ lUJrt = 
£ r [ei cos(ki • r) + ie^ cos(k2 • r)] e~ luJrt , where £ r is the 

lattice field amplitude, €12 = 



2/=F \ \z are field polar- 



izations and ki 5 2 = f 21 (^\[\v =^ yl^) anc ^ Ur = 2nc/X r 
are the the corresponding wave vectors and frequencies, 
respectively. Again following Ref. [43] , one can prove that 
all three (5s4<i) 3 Di(m = 1, 0, —1) states exhibit the same 
ac polarizability a r D (X r ). Based on the data in Refs. [44]- 
[47] . we show a r D (X r ) and the ac polarizability a r P (X r ) 
of the (5s5p) 3 Po state as a function of the wavelength 
A r in the vicinity of 3 fim in Fig. Sl^b). A red-detuned 
magic wavelength X rm ~ 2824.6 nm is located between 
(5s5p) 3 Pi - (bs4d) 3 D 1 and (5s5p) 3 P 2 - (5s4d) 3 Di tran- 
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FIG. SI. Wavelength dependence of the ac polarizabil- 
ities a b ^ r (Xb, r ) (solid lines) and a b p V (Xb, r ) (dashed lines), 
(a) A magic wavelength is located on the blue side of the 
(5s5p) 3 Po — (5s4cf) 3 Di transition, (b) A red-detuned magic 
wavelength is located between (5s5p) 3 Pi — (5s4cf) 3 Di and 
(5s5p) 3 P2 — (5s4cf) 3 Di transitions. In both cases, the ac po- 
larizabilities for all (5s4cf) 3 Di (m = 1,0, —1) states are same. 



and therefore atoms will be trapped in the intensity max- 
ima of the laser field, which permits the straightforward 
construction of a confining dipole trap. 

Note that one can ignore the interference effect be- 
tween blue- and red-detuned laser fields because of their 
extremely large frequency difference. Thus, superimpos- 
ing these two lattice fields creates a setup in which the 
state (5s4d) 3 P and all magnetic sub-levels of (5s4<f) 3 Di 
encounter the same trapping potential, and realizes a ID 
(along the x-axis) long-range interacting many-body sys- 
tem. We would like to point out that this lattice config- 
uration can also be modified to realize a 2D (the x — y 
plane) long-range interacting many-body system. Note 
furthermore, that if only one of the magnetic sub-levels 
of (5s4<i) 3 Di is used experimentally (as discussed in the 
manuscript) one can achieve a magic lattice with less 
stringent laser field geometries. 



sitions, at which is a r P (X r 



ot r D (X rrn ) 



8.3 x 10 3 



a.u., 



[1] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and 

T. Pfau, Phys. Rev. Lett. 9 4, 160401 (2005) 
[2] M. Lu, N. Q. Burdick, S. H. Yo un, and B. L. Lev, |Phys.| 
I Rev. Lett. 107, 190401 (2011)| 
[3J K. Aikawa, A. Frisch, M. Mark, S. Baier, A. Rietzler, 
R. Grim m, and F. Ferlaino, |Phys. Rev. Lett. 108, 210401] 
I (2012)| 

[4] K.-K. Ni, S. Ospelkaus, M. H. G. de Miranda, A. Pe'er, 
B. Neyenhuis, J. J. Zirbel, S. Kotochigova, P. S. Julienne, 
D. S. Jin, and J. Ye, Science 322, 231 (2008) 

[5] M. Saffman, T. G. Walker, and K. M0lmer, 
I Phys. 82, 2313 (2010)] 

|6| H. Weimer an d H~P Biichler, jPhys. Rev. Lett. 105, | 
I 230403 (2010)| 

[7] J. Schachenmayer, I. Lesanovsky, A. Micheli, and A. Da- 
ley, New J. P hys. 12, 103044 (2010). 

[8] I. Lesanovsky, |Phys. Rev. Lett. 106, 025301 (2011)] 



Rev. Mod. 



6 



[9] 
[10] 
[11] 



[121 



[13] 

[14] 
[15] 
[16] 

[17] 



[18] 
[19] 
[20] 
[21] 



[22] 

[23] 
[24] 

[25] 



E. Urban, T. A. Johnson, T. Henage, L. Isenhower, D. D. 
Yavuz, T. G. Walker, and M. Saffman, Nature Phys. 5, 
110 (2009). 

A. Gaetan, Y. Miroshnychenko, T. Wilk, A. Chotia, 
M. Viteau, D. Comparat, P. Pillet, A. Browaeys, and 
P. Grangier, Nature Phys. 5, 115 (2009). 
M. Viteau, M. G. Bason, J. Radogostowicz, N. Malossi, 
D. Ciampini, O . Morsch, and E. Arimondo, Phys. Rev.| 
Lett. 107, 060402 (2011)] 

S. E. Anderso n, K. C. Younge , and G. Raithel, Phys. 
Rev. Lett. 107 , 263001 (2orfj] 

P. Schauss, M. Cheneau, M. Endres, T. Fukuhara, 
S. Hild, A. Omran, T. Pohl, C. Gross, S. Kuhr, and 
I. Bloch, Nature 87, 491 (2012). 

T. Fukuhara, S. Sugawa, M. Sugimoto, S. Taie, and 

Y. Takahashi, Phys. Rev. A 79, 041604 (2009). 

S. Stellmer, B. Pasquiou, R. Grimm, and F. Schreck, 

|Phys. Rev. Lett. 109, 115302 (2012)] 

M. Takamoto, F. L. Hong, R. Higashi, and H. Ka- 

tori, Nature 435, 321 (2005); J. Ye, H. J. Kimble, 

and H. Katori, Science 320, 1734 (2008); T. Akat- 

suka, M. Takamoto, and H. Katori, Nature Phys. 4, 954 

(2008). 

A. V. Gorshkov, M. Hermele, V. Gurarie, C. Xu, P. S. 
Julienne, J. Ye, P. Zoller, E. Demler, M. D. Lukin, and 
A. M. Rey, Nature Phys. 6, 289 (2010); A. V. Gorshkov, 
A. M. Rey, A. J. Daley , M. M. Boyd, J. Ye, P. Zoller 
and M. D. Lukin, |Phys. Rev. Lett. 102, 110503~(2 009) j 
A. J. Daley, M. M. Boyd, J. Y e, and P. Zoller, |Phys.| 
Rev. Lett. 101, 17 0504 (2008) 
R. H. Dicke, Phys. Rev. 93, 99 (1954). 

G. S. Agarwal, Phys. Rev. A 2, 2038 (1970). 
R. H. Lehmberg, Phys. Rev. A 2, 883 (1970). 

H. Zoubi and H. Ritsch, Physica E 42, 416 (2010); J. 
Phys. B 44, 205303 (2011); preprint , arXiv: 1106.4923 
(2011). 



X. Zhou, X. Xu, X . Chen, and J. Chen, |Phys. Rev. A 

81, 012115 (2010)| 

D. F. V. James, |Phys. Rev. A 47, 1336 (1993)| 



[26] 
[27] 



X.-L. Deng, D. Porras, and J. I. Cirac, Phys. Rev. A 77 
03340 3 (2008)| 

P. Hauke, F. M. Cucchietti, A. Muller -Hermes, M.-C. 
Bauls, J. I. Cirac, and M . Lewenstein, |New Journal of| 

Physics 12, 113037 (2010)) 

X. Wang, |Phys. Rev. A 64, 012313 (2001)] 



T. Kennedy, E. H. Lieb, and B. S. Shastry, Phys. Rev. 



Lett. 61, 2582 (1988); K. Kubo and T. Kishi, Phys. Rev 



Lett. 61, 2585 (1988): E. Fradkin, Phys. Rev. Lett. 63 



[28 

[29; 

[30 
[31 

[32; 
[33; 



[34 



[35; 



322 (1989); K. Harada and N. Kawashima, J. Phys. Soc. 
Jpn. 67, 2768 (1998); A. W. S andvik and C. J. Hamer, 
Phys. ^RevTB 60, 6588 (1999)| 

M. Rigol, V. Dunjko, and M. Olshanii, Nature 452, 854 
(2008). 

M. E. Fisher, S.-k. Ma, an d B. G. Nickel, |Phys. Rev.| 
Lett. 29, 917 (1972)1 J. Sak,|Phys. Rev. B 8, 281 (1973)| 
J. M. Kosterlitz, |Phys. Rev. Lett. 37 1577 (1976)| 
M. Picco, preprint , arXiv: 1207. 1018 (2012). 
W. S. Bakr, J. I. Gillen, A. Peng, S. Foiling, and 
M. Greiner, Nature 462, 74 (2009). 

J. F. Sherson, C. Weitenberg, M. Endres, M. Cheneau, 

I. Bloch, and S. Kuhr, Nature 467, 68 (2010). 

F. Pazmandi, G. Zimanyi, and R. Scalettar, Phys. Rev. 



Lett. 75 
S. Das, 



1356 (1995) 
G. S 



Agarwal, and M. O. Scully, Phys. Rev. 
Lett. 101, 153601 (2008)1 
D.-W. Wang, Z.-H. Li, H Zheng 
Rev. A 81, 043819 (2010)] 



and S.-Y. Zhu, Phys. 



[36] T. Tokihiro, Y. Ma nabe, and E. Hanamura, |Phys. RevT 

B 47, 2019 (1993)| 
[37] Note, that in the limit N — >• 00 or in case of periodic 
boundary conditions the two sets of eigenstates coincide 

and thus these narrow lines are abse nt [48]. 

S. D. Jenkins and J. Ruostekoski, |Phys. Rev. A 86, 
031602 (2012)] 

R. Wiegner, J. von Zanthier, and G. S. Agarwal, Phys. 
Rev. A 84, 023805 (2011)] 

B. Olm os and I. Lesanovsky, |Phys. Rev. A 82, 63404 
(2010) 

I. B. Mekhov, C. Maschler, and H. Ritsch, Nature Phys. 
3, 319 (2007). 

I. Lesanovsky, M. van Horssen, M. Gut a, and J. P. Gar- 
rahan, preprint , arXiv:1211.1918 (2012). 
P. Rosenbusch, S. Ghezali, V. A. Dzuba, V. V. Flam- 
baum, K. Beloy, and A. Derevianko, |Phys. Rev. A 79^] 
013404 (2009) 

[44] J. E. Sansonetti and G. Nave, J. Phys. Chem. Ref. Data 

39, 033103 (2010). 
[45] G. Garcia and J. Campos, J. Quant. Spectrosc. Radiat. 

Transfer. 39, 477 (1988). 

J. R. R ubbmark and S. A. Borgstrom, Physica Scripta 
18, 196 (1978)] 
X. Zhou, X. Xu, X . Chen 
81, 012115 (2010)| 
[48] G. Nienhuis and F. Schuller, J. Phys. B 20, 23 (1986) 



and J. Chen, Phys. Rev. A 



